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ABSTRACT 


The wind-accreting X-ray binary pulsar and cyclotron 
line source Vela X-l has been observed extensively dur¬ 
ing INTEGRAL Core Program observations of the Vela 
region in June-July and November-December 2003. In 
the latter set of observations the source showed intense 
flaring - see also Staubert et al. (2004), these proceed¬ 
ings. 

We present early results on time averaged and time re¬ 
solved spectra, of both epochs of observations. A cy¬ 
clotron line feature at ~53 keV is clearly detected in the 
INTEGRAL spectra and its broad shape is resolved in 
SPI spectra. The remaining issues in the calibration of 
the instruments do not allow to resolve the question of 
the disputed line feature at 20-25 keV. 

During the first main flare the average luminosity in¬ 
creases by a factor of ~10, but the spectral shape remains 
very similar, except for a moderate softening. 
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1. INTRODUCTION 


Vela X-l (4U 0900—40) is an eclipsing high mass X-ray 
binary with an orbital period of 8.96437 days (Barziv 
et al. 2001) at a distance of ~2.0kpc (Nagase 1989) con¬ 
sisting of the B0.5Ib supergiant HD 77581 and a neutron 
star. The optical companion has a mass of ~23 Mq and a 
radius of ~30 If . while the neutron star mass is estimated 
to be ~1.8M 0 (Barziv et al. 2001). 

Due to the small separation of the binary system (or¬ 
bital radius: 1.7 R*), the neutron star is deeply embed¬ 
ded in the intense stellar wind (4x 10~ 5 M 0 /yr; Nagase 


et al. 1986) of HD 77581. X-ray line spectra measure¬ 
ments (Sako et al. 1999) show that this wind is inhomo¬ 
geneous with dense clumps embedded in a much thinner, 
highly ionized component. The typical X-ray luminosity 
of Vela X-l is ~4x 10 36 erg/s, but both sudden flux re¬ 
ductions to less than 10 % of its normal value (Inoue et al. 
1984; Lapshov et al. 1992; Kreykenbohm et al. 1999) and 
flaring activity (Kendziorra et al. 1990; Haberl & White 
1990; Kreykenbohm et al. 1999) have been observed in 
the past. 

The neutron star has a spin period of ~283 s (McClin- 
tock et al. 1976). Both spin period and period derivative 
have changed erratically since the first measurement as is 
expected for a wind accreting system. The last measure¬ 
ments with the Burst and Transient Source Experiment 1 
resulted in a period of ~283.5 s. 

The broadband X-ray spectrum of Vela X-l has the typ¬ 
ical shape of accreting pulsar spectra with a power law 
continuum at lower and an exponential cutoff at higher 
energies. This is further modified by strongly vary¬ 
ing absorption which depends on the orbital phase of 
the neutron star (Kreykenbohm et al. 1999; Haberl & 
White 1990), an iron fluorescence line at 6.4 keV, and 
occasionally an iron edge at 7.27 keV (Nagase et al. 
1986). A cyclotron resonant scattering feature (CRSF) 
at ~55 keV was first reported from observations with 
Mir-HEXE (Kendziorra et al. 1992). Makishima et al. 
(1992) and Choi et al. (1996) reported an absorption fea¬ 
ture at ~25keV to 32keV from Ginga. This lower en¬ 
ergy feature has been disputed by BeppoSAX observa¬ 
tions (Orlandini et al. 1998) but supported by phase re¬ 
solved analysis of Mir-HEXE (Kretschmar et al. 1997) 
and RXTE data (Kreykenbohm et al. 2002). 


*See http://www.batse.msfc.nasa.gov/batse/ 
pulsar/data/sources/velaxl.html 
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Figure 1. Sketch of the Vela X-l System during the observations in summer and in winter 2003. The source position along 
its orbit during the INTEGRAL observations is indicated for each INTEGRAL revolution. 



Figure 2. ISGRI image of pointing 013700420010 dur¬ 
ing which the first major Vela X-l flare reached its peak. 
The source was outside of the field of view of the monitor 
instruments. 


2. OBSERVATIONS AND DATA REDUCTION 


As part of the INTEGRAL Core Program the Vela region 
has been observed twice in 2003 for an extended time. 
The first set of observations took place mid June to early 
July (revolutions 81 to 88) the second end of November 
to mid December (revolutions 137 to 141). 

As the observations were not targeted at any specific 
source in order to support different scientific goals, the 
observation strategy was not optimized for Vela X-l. In 
practice this meant that several revolutions were sched¬ 
uled with Vela X-l in eclipse (see Fig. 1) and that for a 
significant fraction of the time the source was not within 
the field-of-view of the two monitor instruments JEM-X 
and OMC. 

The two observation epochs, summer and winter, saw 
a very different behaviour of Vela X-l. While in sum¬ 
mer the source was mostly calm, the winter observa¬ 


tions showed several large flares (Krivonos et al. 2003; 
Staubert et al. 2004), possibly the largest ever observed. 
These flares and the general source behaviour during the 
observations are described in more detail in Staubert et al. 
(2004). Unfortunately the observation strategy explained 
above means that we have no coverage for the flares with 
the monitor instruments. Also the RXTE-ASM did not 
sample this source region during the flares. 

For the summer observations we exctracted long term av¬ 
erage spectra for all three high energy instruments, ex¬ 
cluding times of eclipse and pointings with other prob¬ 
lems like, e.g., radiation belt entries. While these spectra 
only describe an average source state and the exact se¬ 
lection varies from instrument to instrument, the large to¬ 
tal exposure allows to determine a statistically significant 
SPI spectrum beyond 50keV despite the steeply falling 
source spectrum and degradation of the data by a large 
solar flare which occured in revolutions 82/83. The total 
integration times for the combined spectra are 340.2 and 
333.2 ksec for JEM-X and ISGRI. respectively, where 
concurrent observations were chosen and 775.6 ksec for 
SPI. 

For the winter observations we have concentrated on the 
spectra before and during the first big flare. In this case 
we have generated pulse phase resolved spectra covering 
the rise, the peak and falling flank and the following off- 
pulse region connected with the so-called “main pulse” 
as defined, e.g., by Kreykenbohm et al. (2002). The used 
bins are indicated in Fig. 3. 

The data reduction was done with the OSA3 release 
of the scientific analysis software and the correspond¬ 
ing calibration files. Due to the remaining problems in 
the ISGRI response matrix with “snake-like” residuals, 
we assumed systematic uncertainties of 10% for the ex¬ 
tracted ISGRI spectra. For JEM-X a known problem ex¬ 
ists in reconstructing fluxes below ~ 10 kcV, in our case 
this leads to artificially enhanced intrinsic absorption in 
the, typically strongly absorbed, source. 
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Figure 3. The phase intervals chosen for the pulse phase 
resolved spectroscopy. 


Since no direct support for pulse phase resolved anal¬ 
ysis is yet available in the analysis software, we con¬ 
structed specific GTI information for the different phase 
bins. These GTIs are based on the barycenter and bi¬ 
nary corrected pulse period determination (Staubert et al. 
2004) from which then times in the instrumental time 
frame were derived. 


3. RESULTS 


The combined spectra of JEM-X and SP1 for the sum¬ 
mer observations can be well fitted at lower energies by a 
power law model with exponential cutoff (XSPEC model 
CUTOFFPL), if one includes free relative normalization 
between the two spectra, to allow for (inter-)calibration 
uncertainties. A powerlaw with a Fermi-Dirac cutoff 
(Tanaka 1986) works equally well and gives very similar 
results, but the onset of the cutoff is not well determined, 
thus we use the simpler model here. 

Above 50 keV there is an evident structure in the fit resid¬ 
uals for SPL see Fig. 4. This structure can be fitted by in¬ 
cluding a cyclotron scattering feature at Efe54keV with 
rr«7 keV. Since the SPI energy resolution at these ener¬ 
gies is ~1.6keV (Attie et al. 2003), the line shape can 
actually be resolved as far as statistics permit it. The line 
parameters do not depend siginificantly on the exact con¬ 
tinuum - simple exponential cutoff or Fermi-Dirac cutoff 
- used. 


Combining JEM-X and ISGRI spectra, the missing flux 
above —50 keV is again clearly visible in the fit residuals 
- see Fig. 5 - and can be fitted with a line at i?«52 keV. 
Combining spectra from all three high energy instruments 
does not allow a good fit due to calibration differences 
between ISGRI and SPI even though the fit parameters 
are mostly compatible (see Tab. 1). 

Unfortunately, the current calibration and cross-cali¬ 
bration uncertainties for the INTEGRAL instruments do 
not allow any firm statement about the existence or not 
of a —25 keV line in the spectra. All spectral modeling 
has been done without including such a feature but a fea¬ 
ture of the strength reported, e.g., by Rreykenbohm et al. 
( 2002 ) would be fully consistent with the data at the mo¬ 
ment. 
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Figure 4. Combined JEM-X and SPI spectra averaged 
over the good pointings during the summer observations. 
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Figure 5. Combined JEM-X and ISGRI spectra averaged 
over the good pointings during the summer observations. 


Table 1. Selected best fit parameters for the combined, 
phase averaged JEM-X & SPI spectra and JEM-X & 
ISGRI spectra respectively using the XSPEC CUTOFFPL 
model. Note that the values are still preliminary due to 
calibration uncertainties. 


Parameter 

JEM-X + SPI 

JEM-X + ISGRI 

photon index 

0 +°- 10 
u - ja -0.20 

0 43 + 0.10 

folding energy 

11.3i?;^keV 

12 . 1 + 0 ^ keV 

line center 

53.6±f;| keV 

51-9±|i keV 

line a 

7-3±i;f keV 

4 . 0^ ; 5 keV 

line depth 

0 63 +0 13 
U - OJ -0.08 

0-38l^ 16 

Xred with line 

1.19 

1.02 

Xred w/o line 

1.85 

1.17 



























ISGRI: main peak rise 


ISGRI: main pulse peak+fall 


ISGRI: off pulse 
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Figure 6. Comparison of ISGRI spectra before and during the first large flare of Vela X-l. From left to right the panels 
show the spectra for the rising flank of the main pulse, the peak and falling flank and the following off-pulse region. In 
each panel the lower spectrum is that of the pre-flare while the upper is taken over the flare. The spectral softening is 
directly visible. 

Table 2. Selected best fit parameters for phase resolved ISGRI spectra before and during the large flare of Vela X-l in 
revolution 137, using the XSPEC CUTOFFPL model. Note that the best fit values are still preliminary. 

Main pulse rise Main pulse peak & fall Off pulse 


Parameter 

pre-flare 

flare 

pre-flare 

flare 

pre-flare 

flare 

norm 

Q 1 Q + 0.01 
u - A :> - 0.01 

2 14+ 0 - 62 
z * ah - 0.05 

0 21 + 0 - 03 
w ‘ Zi - 0.02 

r\ QQ + 0.09 

0 09+ 0 - 03 
u - uy -o.oi 

1 4^+ 0 - 25 

A - HO — 0.29 

folding energy 

11.2+1 keV 

9.6iS‘i keV 

10.0±g;| keV 

8-9±o.i keV 

9.6+;? keV 

8.5+;® keV 


SPI: main pulse rise 


SPI: main pulse peak+fall 


SPI: off pulse 





Figure 7. Comparison of SPI spectra before and during the first large flare of Vela X-l. The panels follow the same pattern 
as in Fig. 6. Due to the relatively short integration times, the feature at ~50 keV is not clearly visible. 

Table 3. Selected best fit parameters for phase resolved SPI spectra before and during the large flare of Vela X-l in 
revolution 137, using the XSPEC CUTOFFPL model. Note that the best fit values are still preliminary. 

Main pulse rise Main pulse peak & fall Off pulse 


Parameter 

pre-flare 

flare 

pre-flare 

flare 

pre-flare 

flare 

norm 

0-25+;°® 

9 17 + 0-81 

Z * A '- 0.58 

0-321+05 

4 16~*~ 2 ' 50 

^• AO — 1.56 

n i 3+ 0 - 09 
u - a :) - 0.05 

1 ^ Q + 2 - 10 

a - j ^- 0.90 

folding energy 

11 -4±g;I keV 

10.2±^;J keV 

11-5+8 keV 

8.3+o keV 

11.0+8 keV 

9.6+;? keV 




























































































3.1. Winter observations 

For the pre-flare and flare data of the winter observations 
we are restricted to spectra from the two main instruments 
as explained in Section 2. Fig. 6 and Fig. 7 show for com¬ 
parison the phase resolved spectra of before and during 
the first big flare. While the flux increases by a factor of 
about 10 during the flare, the spectral shape stays more 
or less the same, except for a certain softening, which is 
directly visible in the count spectra. 

Using the same model as for the summer observations the 
folding energy parameter, which can be seen as a measure 
of temperature, decreases by 1-1.5 keV during the flare, 
see Tab. 2 and Tab. 3. For these fits the photon index 
of the continuum and the line width and energy were ill 
determined and fixed to the values found for the summer 
data, see Tab. 1. This had little effect on the continuum 
parameters. 


4. DISCUSSION 


While the remaining uncertainties in the calibration of the 
individual instruments as well as in the intercalibration do 
not allow any firm statement about a possible line feature 
at 20-25 keV we clearly observe and resolve - thanks to 
the high spectral resolution of SPI- the well-known line 
feature above 50keV. The half-width of 7.2^2 keV ob¬ 
tained from SPI data is in the typical range of line widths 
reported from observations with other instruments and 
expected from thermal broadening. It falls inbetween the 
relatively narrow line solution reported by Coburn et al. 
(2002) and the very broad line width given in Orlandini 
et al. (1998). It should also be kept in mind that this line 
width is obtained for a long-term time averaged spectrum, 
as for shorter integration times the statistics do not allow 
to constrain the fluxes in this energy range anymore. 

Based on the theoretical cyclotron line shapes proposed, 
e.g., by Araya & Harding (1999), the broad, flat shape 
and lack of additional structure could be taken as indica¬ 
tion that this line is indeed a harmonic and not the more 
structured fundamental. Alternatively, this would indi¬ 
cate observation at a relatively narrow angle to the mag¬ 
netic field and preferably of a cylindrical emission region. 

It is interesting to see that despite a rapid, strong brighten¬ 
ing, the source spectrum appears to remain very similar 
to that obtained during more ususal flux levels. With a 
normal luminosity of ~4x 10 36 ergs/s the wind-accreting 
source remains in the sub-Eddington regime of accretion 
even during the flares. Given the well known pulse to 
pulse variations (Staubert et al. 1980) one can speculate 
that the normal mode of accretion for Vela X-l is a rel¬ 
atively sparse and inhomogenous flow which does not 
fill the accretion column completely. In this picture the 
spectral softening observed could be caused by additional 
comptonisation of X-ray photons by a more densely filled 
accretion column. 

A more detailed analysis of this extensive and higly in¬ 
teresting dataset than can be done for these proceedings 


is surely called for in the future. With further improved 
flux reconstruction and calibration information not only 
the question of the 20-25 keV line can be settled but we 
will also arrive to do detailed quantitative studies of the 
spectral evolution in calm and flaring times. 
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